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Abstract: We provide a dual gravity description of heavy atomic nuclei, via 
AdS/CFT correspondence. In holographic QCD such as Sakai-Sugimoto model, 
baryons are D-branes wrapping a sphere in 10 dimensional curved spacetime, so 
any nucleus is a collection of A such D-branes where A is mass number of the nu- 
cleus. Quantum theory on the nucleus is ADHM-like U{A) Yang-Mills-Higgs theory 
on the sphere. Taking a large A limit (corresponding to heavy nuclei) leads to a 
dual gravity describing collective excitataions of constituent nucleons of the heavy 
nucleus. This dual gravity computes spectra of the heavy nucleus, and gives discrete 
states whose gap roughly agrees with experimental nuclear data. 
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1. From superstring to nuclear theory 

Application of a superstring technique, the gauge/gravity (AdS/CFT) correspon- 
dence in, 0, to low energy QCD has provided quite remarkable progress on hadron 
physics. This subject called holographic QCD has grown up to be a major research 
arena of string theory. Dual gravity description has revealed various aspects of low 
energy QCD which were unreachable by conventional analytic methods because of no- 
torious strong coupling. They include spectra of glueballs, mesons and baryons, and 
interactions among them, and even phase transitions at finite temperature/density. 
However, almost all of the results are on hadron physics,* not really nuclear physics, 
dare to mention. In this paper, we take a first step toward nuclear physics: we 
provide a dual gravity description of heavy nuclei. 

The essence of the gauge/gravity correspondence is the large N limit, where N 
refers to the rank of the gauge group U{N) of the gauge theory living on N coincident 
D-branes. This limit, together with large A ('tHooft coupling) limit, allows a dual, 
holographic, equivalent gravity theory on a near horizon geometry of black brane 
solutions created by the N D-branes. 

Baryons in QCD-like gauge theories are described, in the dual gravity description, 
by in fact additionally put D-branes wrapping a higher dimensional sphere in 10 
dimensional curved spacetime [§]. Any nucleus is a collection of baryons (nucleons), 
so in holographic QCD the nucleus is a collection of D-branes. Therefore, heavy 
nuclei with large mass number A can have a dual gravity description provided by 

* Inclusion of baryon chemical potential or finite baryon density is a hot topic in holographic 
QCD, see for a partial list of the references. 
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a near horizon geometry of those "baryonic" A D-branes, in large A hmit. In this 
paper we pursue and reahze this idea. We describe a possible dual geometry, check 
its validity as a limit of string theory, and provide a comparison with nuclear data. 

In the next section, we describe a gauge theory on the A D-branes representing 
the nucleus, with a brief review of baryons in AdS/CFT. Then in Sec. 3, we take 
a large A limit and give a gravity dual, whose fluctuation analysis is discussed in 
Sec. 4 to provide spectra of heavy nuclei. The last section is for a conclusion and 
discussions on possible relations to shell model and liquid drop model of nuclei, and 
black holes. 



2. Gauge theory on nucleus 

Among various holographic models of QCD, Sakai-Sugimoto (SS) model P] is a 
successful one given by top-down approach from string theory. We will work in this 
SS model. Baryons were introduced in the model in @, 0, while in the context of 
AdS / CFT correspondence description of baryons in terms of D-branes has a much 
longer history 

For A/" = 4 U{N) supersymmetric Yang-Mills theory, the gravity dual is type IIB 
supergravity theory on AdS^ x S^. Baryon is described by a D5-brane wrapping the 
S^. All the spatial worldvolume dimensions of the D5-brane are in this "internal" S^, 
thus it is point-like in the AdS^. The reason why this D5-brane is a baryon is that 
this supergravity has a background 5- form field strength penetrating the S^, with 
units of flux. This induces electric charge A^ on the spherical D5-brane. Consistency 
of the charge conservation requires fundamental strings emanating from the D5- 
brane surface since electric charges are created at the end points of fundamental 
strings on a D-brane. The strings ending on the D5-brane should be elongated to 
the AdS boundary, and those are identified as (infinitely heavy) quarks. Thus the 
spherical D5-brane consists of A^ quarks, that is, a baryon. Precise supersymmetric 
configuration of the D5-brane with the electric charge was given in [Q. 

One can collect A D5-branes on top of each other. On the worldvolume (~ 
R X S*^) of the A D5-branes, as usual for any D-brane systems in string theory, there 
appears a U{A) super Yang-Mills theory. It includes scalar fields in adjoint represen- 
tation which measure relative position of the D5-branes. Once the worldvolume 
is decomposed into spherical harmonics, the lowest mode depends effectively only on 
time, so the theory reduces to a f/ (A) supersymmetric matrix quantum mechanics. 

The situation is qualitatively similar for the SS model. The model provides 
various interesting features of baryons [|^, 0, ^ As for the gravity dual of the 
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gluonic part of QCD, the model adopts so-called Witten's geometry |[Tl[| , 

ds' = (^y\r^dx^dx'^ + f{U)dr') + I^^Y' ^ (2.1) 



'f;\3/4 



(2.2) 



with a 4-form field strength fiux penetrating the S*^ specified by the volume element 
dQ^ above. This is a supersymmetry-breaking solution of type IIA supergravity in 
10 dimensions. The supergravity parameters are related to QCD parameters as 



R 



V 2Mkk 



1/3 



Ukk — gS'YM^c^KK^s, ds 



9ym 



2ttMkkIs 



[2.3) 



Flavors are introduced as Nf D8-branes in this geometry. The D8-branes extend 
all spatial directions except for r parameterizing 5*^ whose radius is 1/Mkk- All 
the hadronic physical observables computed in this model are given as functions 
of Nc,gYM and Mkk- Because the geometry ( pTTP is truncated smoothly at ?7 = 
Ukk, one can introduce a new coordinate {z,y) instead of {t,U), a.s z + iy = 
\J (U^/Ukk) — U^^exp{iMKKT), then the D8-brane location is expressed as y = 0. 
The D8-brane gauge group U (Nf) corresponds to unbroken global symmetry in QCD, 
the vector part of the chiral symmetry. 

The baryons in the SS model are given by D4-branes wrapping the 0, in the 
gravity dual side. See the table below, and Fig. |l[ It is particle-like and localized in 
the {x^,x',x^) space. 
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The D4-branes sitting on the flavor D8-branes can be represented by Yang-Mills in- 



stantons (localized in (x^, x^, x^, z) directions) of the D8-brane gauge theory [Q, and 



in fact, due to a Chern-Simons coupling on the D8-branes, the Yang-Mills instanton 
sources an overall U{1) gauge field on the D8-branes and thus possesses a baryon 
number. 

Collecting the baryon D4-branes, we consider A of them to describe a nucleus. 
On the A D4-branes wrapping the S*^, one has 1+4 dimensional U{A) gauge theory 
on R X 5*^. At low energy of these D4-branes, we have U{A) gauge fields, adjoint 
scalar fields $i {i = 1, 2, 3, y, z) and fundamental scalar fields (a = 1, 2, ■ ■ ■ , Nf). 
The adjoint and fundamental scalar fields are responsible for the location and the 
size of the instantons, respectively. The theory is a deformed version of "ADHM" 
gauge theory: it is well-known that the theory on BPS Dp-branes in the presence 
of D(p + 4)-branes gives rise to ADHM equations of Yang-Mills instantons as BPS 
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A D4s 



equations |]T3|, [T^. In our case, the theory is not supersymmetric, due to the su- 
persymmetry breaking in the bulk curved geometry. Therefore the theory on the 
D4-branes has nontrivial potential terms which are difficult to analyze. We consider 
generic mass number A for a nucleus. In particular when A is large, the theory is 
complicated. 

The case of a single baryon was an- 
alyzed in 0. In fact, there appear two 
deformation terms; one is due to the 
background curved geometry which gives 
rise to a potential for the $2 field, and 
the other is due to the background 4- 
form flux providing the Chern-Simons 
coupling. The size of the instanton pa 
obtains a potential from both of these 
two effects. Larger $2 or larger pa cost 
energy due to the background metric. 
So $2 is stabilized at 2; = 0, i.e. the 
instanton localizes at the tip of the ge- 
ometry. When the size pa is small, elec- 
tric charge (=instanton number) density localizes and costs more energy, so the size 
is stabilized at a certain finite value. This was computed in H as 




f/KK U 

Figure 1: Baryons are D4-branes sitting at 
the tip of the geometry (depicted by a shaded 
blob). They wrap which is not shown in 
this figure. 



p = 3v37r 



1/4 



(2.4) 



KK 



Quantization of this electrically charged (=dyonic) instanton on the D8-branes was 
carried out in to obtain a spectrum of baryons. 

Normally, in holographic QCD, all the KK modes on S*^ are discarded, simply 
because those excitation don't exist in QCD. If we apply the same philosophy to 
the U{A) gauge theory on the baryon D4-branes wrapping the S"', the gauge theory 
on the D4-branes reduce to a matrix quantum mechanics. This theory is, in the 
holographic QCD, a microscopic theory of multi-body nucleons. 



3. Large A limit and gravity dual 

Since the U{A) gauge theory is difficult to treat, we take large A limit and consider 
its dual geometry. 

Let us explain the whole picture first. Note that the U{A) gauge theory lives in 
the gravity dual of QCD, already. The A D4-branes are in the curved background 
created by the Nc color D4-branes. In this situation, we consider a large A limit. 
The large number of D4-branes in the Witten's geometry back-reacts and deforms the 
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geometry. We take the near-horizon hmit of this back-reacted geometry by getting 
closer to the baryon D4-branes. So, effectively, we take the AdS/CFT duality twice. ''^ 

Before taking the near horizon limit, we need a supergravity solution. In fact, 
it is not just a supergravity but the one with the flavor D8-branes. The system 
consists of the supergravity in Witten's background + U{Nf) gauge theory on the 
D8-branes. We introduce A dyonic instantons on the D8-brane gauge theory in the 
Witten's background, and need to solve the equations of motion for the gravitating 
dyonic multi-instantons there. 

This turns out to be very difficult, in various sense. First, multi-instanton solu- 
tions with large instanton number are difficult to construct (for example, see ||15|| ). 
Second, the instanton equation has two deformations described above. Thirdly, the 
instantons live only on the D8-branes which are localized in the bulk, so the gravity 
problem is essentially similar to that of braneworld black holes, which is known to 
be difficult to solve even numerically. 

So, we need to assume some gravitating dyonic instanton 
solution, to proceed. For that, let us study expected property 
of the solution. After all, we come to a conclusion that the 
solution is approximated by asymptotically-flat black-4-branes 
whose singularities are distributed on a 3-ball with a certain 
radius. See Fig. |^. 

First, let us recall a well-known property of nuclei. Among 
nucleons a very strong repulsive force exists, while there is a 
binding force due to light meson exchange. Due to this, basi- 
cally nucleons are stabilized at a certain distance. Resultantly, 
any nucleus has a common nucleon density, described by a 
nucleus radius formula 




Figure 2: The 3-ball 
distribution of A D4- 
branes (shaded re- 
gion) on D8-branes. 



Rn{A) = pA 



1/3 



(3.1) 



Experimentally this is valid well, with p ~ 1.1 [fm]. Interest- 
ingly, the value ( p.4[ ) obtained in 0, with Mkk ~ 500 [MeV] which fits the baryon 
spectrum, coincides with the experimental value of p. It should be possible to repro- 
duce this nuclear force from the SS model |T^, but here we assume this homogeneous 



distribution of dyonic instantons over a 3-ball on the D8-brane. The radius of the 



^This is not rare situation. For example, Coulomb branch oi Af — 4 super Yang-Mills was 
studied by a probe D3-brane in AdS^, and if one collects many D3-branes there, one needs to go to 
a back-reacted geometry which is a multi-center D3-brane solution . What we are doing here is 
qualitatively similar to that. 
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3-ball is, due to 

^1/3 

Rp, ~ . (3.2) 

We neglect all numerical coefficients in this paper. Note that the ball is in {x^,x'^, x^) 
space, and the dyonic instantons are localized at z = 0. These 4 directions are spatial 
directions of the D8-branes except for the directions. The dyonic instantons feel 
a potential along z and are stabilized at z = 0. In fact dyonic instantons with 
wave functions which expands along z direction correspond to excited nucleons 
which are not expected to be present in heavy nuclei. 

Second, notice that the size of the nucleus R^ ( p.2|) is small for large A. We 
are interested in the near horizon limit of the A D4-brane geometry, so only a small 
region of the bulk geometry is necessary. The D4-branes are at 2; = 0, so, before the 
back-reaction, the geometry is approximated by 

ds' = (^^^ '^^^''dx^dx'' + ^ '^\dy' + dz') + R'/^Ullidnl (3.3) 

This is almost fiat space, with rescaled variables x'^ = {Ukk/ RY^^x'^ and {y,z) = 
{2/?,){R/UKKf/^{y,z). We place A D4-branes distributed in a 3-ball with its radius 
Rn (|3.2|). If we neglect the presence of the D8-branes, the curvature of the 5*^ and 
also the background flux, the back-reacted solution is just a multi-center BPS black 
D4-brane, 

ds^ = r''H-{dx^f+ lim R%dnl\ + f'^{dx'f + df + d-z") , (3.4) 

\ ijg4— >oo / 

+ rj^\T-wjl (3.5) 

7=1 

with r = {x\y,z) {i = 1,2,3), and wi are transverse location of the D4-branes 
(labeled by / = I,-- - ,A). Here g'^ is the string coupling constant evaluated at 
U = Ukk (equivalently, z = y = 0), 

When the D4-brane distribution forms a uniform 3-ball of radius R^, 



atiR% Jb^ n 



/ = 1 + / d'w (3.7) 



We need to take a near horizon limit. For this, it is necessary for R^ to be 
smaller than "would-be-AdS radius" Rq which is determined by the factor in the 
harmonic function / in ( p.5|) , 

Rq = {nA g'j'/%, Rq = Ro/y^i- (3.8) 

we use charge radius of proton (which was obtained in the SS model in instead of ( |2.4| ), 
we have a different scahng in A. The author thanks S. Sugimoto for pointing this out. 
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If the condition is violated, the ball is larger than the "AdS radius" and thus 
AdS/CFT correspondence doesn't work. In (|3.8|) , we have included the rescaling 
factor ^Jgii = {Ukk/ RY^^ necessary for scaling x back to x. Substituting (|2.3| ), we 
obtain 

'? 41/3 

^0 = ^^7^ • (3.9) 

The condition Rjq < Rq is satisfied if 

Nl'"^ <V\. (3.10) 

This can be satisfied for a chosen order of large Nc and large A limit. ^ (If we substitute 
values used in which are fixed by fitting the experimental p meson mass and pion 
decay constant, this inequality is satisfied.) In this parameter region, the near horizon 
limit makes sense. 

4. AdS/CFT dictionary and nuclear spectrum 

According to the standard AdS/CFT dictionary, gauge invariant operators of the 
U{A) gauge theory on the nucleus correspond to fluctuations of supergravity fields 
and fields on the flavor D8-branes, in the background described in the previous 
section. More precisely, for example, the dictionary consists of 

tr([$i,$fc][$j,$fc]) + --- ^ Sgij 

phipb + ■■■ ^ 5iAi)\ 

Fluctuation spectrum of the bulk fields is the spectrum of these composite operators. 
In particular, collective motion of the constituent nucleons is expressed by these 
bulk fields. For example, the collective motion involving movement of the nucleons 
corresponds to gravity fiuctuations, while the gauge field fiuctuation on the D8- 
branes are mainly responsible for nucleon density waves in the nucleus. Bulk dilaton 
fluctuation includes all the effects. 

We can perform a supergravity fluctuation analysis using the explicit metric 
given in the previous section. The simplest choice of the fluctuation fleld is the 
dilaton fluctuation, which satisfles free equation of motion in the curved geometry,* 

^ -.dMe-"t'y^g''''dr,<P = 0. (4.1) 
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§It was pointed out in JTzt (and argued also in |^) that the original supergravity approximation 
for the background (2.1) is valid for ^ g'^ ^jv • This contradicts (3.10). However, note that 



these conditions are just for order estimates which neglect numerical factors. In this paper we 
simply assume that there is a parameter region (for gyMdndNc) at which both the above condition 
and the condition Rpf < Rq can be satisfied. The author thanks H. Ooguri and S. Sugimoto for 
bringing this discussion to him. 

*Generically fluctuation of dilaton can mix with other gravity fluctuations. Here we are interested 
in only the order of gap of discrete energy spectrum, as we will see below. 
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It turns out that, to solve this equation in the background black 4-brane solution, 
one needs to do a numerical calculation, even if the presence of the D8-branes is 
neglected.^ 

What we want to see is if the spectrum is discrete or not, and if it is, what is the 
approximate gaps appearing in the spectrum. For this purpose, we shall look at a 
similar black-brane system for which extensive analyses have been performed already 
in the past. One very similar black brane solution is found in the supersymmetric 
context, that is. A/" = 4 supersymmetric Yang-Mills theory in Coulomb branch. 



In [|T9|, |2^, ^ (see also [[16| , |2^), certain ball-like distributions of D3-branes are 
considered. It was found there that dilaton fluctuation develops a mass gap of order 
Rj\f/Rl where in that case R^, Rq are ball radius and the AdS radius, respectively. 
In particular, it is necessary for the ball dimension to be larger than 2, in order to 
have a discrete spectrum |T^. The typical scalar fluctuation spectrum |]19| is-'- 



E^r^^n n = l,2,--- (4.2) 



The striking point is that the continuous distribution of the D-branes makes the 
spectrum discrete and gives a mass gap. The result is robust qualitatively, for 3- and 
5-balls, and 3-sphere.^ 

Interestingly, ( [4.2| ) has apparent similarity to the case of gravity background for 
confining phase. For confining gauge theories, geometry is truncated at a certain 
radius and one cannot go deeper beyond the wall, which gives rise to the discrete 
spectrum. In the present case, this wall is the ball surface. 

Our background is made by the A D4-branes while the above is for distributed 



D3-branes. However, there is an argument |^ for qualitative understanding of 
the spectrum formula ( [4.2| ). The denominator of ( [4.2| ) is i?7v which is naturally 
understood as energy coming from a string stretched between the distributed D- 
branes. The average length is proportional to Rn- In order to fix the dimensionality, 
one needs l/R^ in the formula. 

Assuming this robustness of the formula ( |4.2| ) for distributed D-branes in a ball, 
we can apply it to our particular case. If we substitute ( p.2| ) and (|3.9|), we obtain 



E^r-^j^j^MKKn (n=l,2,---) (4.3) 



^Precise treatmnt is under investigation [ p^ . 

tin this paper -we consider Sakai-Sugimoto model, but we expect that starting from other models 
may lead to similar results. For D3-D7 model (or for cut-ofF AdS^ -with hard -wall in bottom-up 
approach |2j]), baryons are D5-branes -wrapping S^, so the heavy nucleus is dual to near horizon 
geometry of A D5-branes. 

§For a 2-ball (a disk), there is a mass gap of order Rjq / Rq -while the spectrum is continuous 
above it 1191. 
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We neglected all the numerical factors. This (|4.3| ) is the excitation spectrum of heavy 
nuclei with mass number A, as a result of AdS/CFT correspondence. 

It is surprising that this formula is actually consistent with the experimental 
data of heavy nuclei, as follows. It is known from experiments that heavy nuclei have 
coherent excitations (phonons) of constituent nucleons, called giant resonances,^ 

En = uj{A)n (n = l,2,---) (4.4) 

Excitations with lower n have been analysed in detail experimentally. Phenomeno- 
logical models, including liquid drop model of heavy nuclei, indicates consistent har- 
monic behaviour in n. 

Interestingly, the harmonic-oscillator-like spectrum ( [4.4| ) proportional to n is sim- 
ilar to what we obtained, (^73|) , via gauge/gravity duality. In the liquid drop model, 
the harmonic spectrum emerges from nuclear surface oscillation, so the spectrum is 
considered to arise from collective motion of constituent nucleons. Our AdS/CFT 
interpretation, due to the dictionary, is consistent with this. 

Finally, let's compare numerical values of our formula with experiments. We are 
interested only in order of magnitudes. The gap among discrete states in our formula 
is evaluated as 



Rn Nc 



2/3 

Mkk ~ 0(10^) x A^^/^ [MeV] (4.5) 



Rl A^/^Vx 

with typical values of the parameters in the SS model, Mkk ~ 500 [MeV], A ~ 
0(10) — (9(100) and = 3. Experimental results for giant resonance of heavy 
nuclei with A > 60 shows an empirical formula, E ~ 80 A~^^^ [MeV], for the first gap 
of O"*" isoscalar excitation. We considered a dilaton fluctuation which corresponds to 
O"*", and our result ( |4.5|) is not so much different from the nuclear data. In particular, 
we reproduced the A dependence of the resonance excitations. As we adopted very 
crude approximation (we discard all numerical factors and even adopted similar but 
different geometry), one shouldn't take the numerical values seriously. 



5. Conclusion and discussion 

In this paper, we have given a dual gravity description of heavy atomic nuclei, by 
applying the gauge/string duality (AdS/CFT correspondence). Since nucleons are 
described by D-branes wrapping a sphere in curved geometry of holographic QCD, 
on a nucleus with mass number A there appears a U{A) gauge theory. We took a 
large mass number limit A ^ oo. Dual gravity description is valid in this limit, and 

^The author is indebted to T. Nakatsukasa for invaluable discussions on this correspondence, 
and also to T. Matsui for his helpul explanations. 
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we obtained a near horizon geometry corresponding to the heavy nucleus. The cor- 
responding supergravity solution has discrete fluctuation spectra, and we compared 
them with nuclear experimental data. 

Some discussions are in order. We made some arguments for approximating the 
dual geometry, and it is important to validate those by flnding an explicit geometry 
and performing fluctuation analysis. More precise treatment may reveal other aspects 
of the fluctuations. For example, we can put a probe D4-brane as an additional 
nucleon, and consider a probe dynamics in this near horizon geometry of A D4- 
branes. This may lead to a shell model potential of nuclear theory. 



Since our result ([4.3|) has an A dependence, one can compare it with A depen- 
dence in experiments. In the standard holographic QCD, N of "large N" Maldacena 
limit is flxed to be A^^c = 3 for comparison with nature. On the other hand, we 
have another number A counting the nucleon D-branes, which appears in nature 
as a variety of heavy nuclei. This enlarges possibility of comparison of AdS/CFT 
with nature. Another aspect which can be distinguished from other applications of 
AdS/CFT to nature is that clearly we have a U{A) gauge theory on the CFT side. 
We don't rely on any universality for conformal physics. 

Since the supersymmetries are completely broken, one may wonder if event hori- 
zon may develop in the dual geometry and the dual gravity would become a black 
hole. We don't expect it, because, if this were the case, the fluctuation spectrum 
would become continuous (as known in deconflnement phase of holographic QCD). 
However, since black holes are described by fluid dynamics holographically ||25|, one 
can speculate that the liquid drop model of heavy nuclei may be related to dual 
geometries through the holographic hydrodynamics. In fact, dissipation of excita- 
tions on a nucleus is a target of research for many decades." Investigation along this 
direction is in progress. 
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